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Abstract. Electrical power system calculations rely heavily on the Y^us matrix, which 
is the Laplacian matrix of the network under study, weighted by the complex-valued ad¬ 
mittance of each branch. It is often useful to partition the Ybus into four submatrices, to 
separately quantify the connectivity between and among the load and generation nodes 
in the network. Simple manipulation of these submatrices gives the F^q matrix, which 
offers useful insights on how voltage deviations propagate through a power system and on 
how energy losses may be minimized. Various authors have observed that in practice the 
elements of Elg are real-valued and its rows sum close to one: the present paper explains 
and proves these properties. 
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1 Introduction 


Currents (!) and voltages (U) in an electrical power system are related by the admittance 
matrix, , which is generally constructed to have the properties of a weighted Laplacian 
matrix (disregarding shunt element modelling), see [T], [5]. It can usefully be partitioned 
as follows by reordering to group generator (G) and load (L) nodes separately: 


Ig 

II 



( 1 ) 


Where Ig,Vg G , Il,Vl & C" and 


Ygg Ygl 
Ylg Yll 


Where Ygg G Yll G and Ylg = Yq^ g C”^™. With ()^ we denote the 

(non-conjugate) transposed tensor. The Y^us matrix is square by definition. It has N 
rows and N columns, where N is the number of nodes in the power system. The system 
is reordered and partitioned to group G buses and L buses separately. The n L buses 
plus the m G buses equals the total numbers of buses N. Although there may be buses 
in a power system that don’t connect any L or G but are just passive interconnection 
points, in this analysis they would be grouped with the load buses. Typically Yll will be 
larger. For instance, on the IEEE 118 bus system we have 19 G nodes and 99 L nodes, 
i.e. Ygg G and Yll G The other sub-matrices are not generally square. 
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in this case Y^c G and Yql S Note that Y^g = ^gl- The justification 

for this assertion is the guaranteed symmetry of the Yb„s matrix (it is constructed to 
have this property, which is only broken under the rare circumstance that phase-shifting 
transformers are included in the model). The Y^us matrix is complex-valued. The physical 
properties of electrical conductors imply that the imaginary part of each element is larger 
than the real part in nearly all cases. In the below figure is plotted the real and imaginary 
component of every element of the Y^^g matrix for the 118 bus system. The imaginary 
component is on the vertical axis: The real:imaginary ratio is approximately equal for all 
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Figure 1: Empirical Real-Imaginary ratio 


elements, as can be seen in Figure 1. 

The partitioning of the Yb„s matrix into load and generation blocks was introduced by 
Kessel and Glavitsch in [3], and has subsequently been applied to a wide range of power 
engineering problems g], 0, [S], [7], [5], 0, gBi P^ - Recent work by Sikiru et. 
al. [13], [U, im has used the partitioned Y^us approach to develop a more fundamental 
understanding of a power systems inherent connective structure, using, for instance, Schur 
complements and eigen analysis. Relatedly, recent work by Abdelkader et al. [16) . m, 
[18] offers a deeper conceptual understanding of these Yi,us partitions, demonstrating how 
they allow the separation of currents in the network into load and generator induced 
components. Notably [18] demonstrates how a strictly equal real:imaginary ratio for every 
element in the matrix brings one component of the network physical power losses to 
zero. 

From Q: 

Ig = YqqVq + YglVl ( 2 ) 

and 

II = YlgVg + YllVl- (3) 

Rearranging 

Vl = ZllIl — ZllYlgVg- (4) 
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Where 


\ if det{YLL) = 0 / ■ 

The matrix Y^^ is the Moore-Penrose Pseudoinverse of Y^Lj calculated by the singular 
value decomposition of Yll, see US], HO]. Substituting for Vl in 


Ig = O^GG — YglZllYlg)Vg + YglZllIl- (5) 


Equations Q and ([^ are typically represented in matrix form, which permits useful 
engineering applications: 


Vl ' 


Zll 

Flg 

■ II ■ 

Ig 


Kgl 

Yggm 

Vg 


Where: 

Yggm = Ygg — YglZllYlg 


and 

Flg = —ZllYlg = —Kgl ( 6 ) 

Various works (e.g [T4] . a, ini) have noted, and, indeed, relied upon, the observation 
that the elements of are in practice real-valued and that its rows sum to unity, or 
close to. 


2 Main results 

As written in the previous section, the ratio (/to : Re) of each entry of the Y^us matrix will 
in practice tend to be fairly homogeneous though this isn’t always guaranteed. The Yf,us 
matrix will have generally diagonal elements that are positive real and negative imaginary. 
Off-diagonals will be negative real and positive imaginary. Based on these observations 
we can provide the following Proposition. 

Proposition 2.1. Assume the Yb^s matrix as defined in Q. If the entries in each 
row have the reahimaginary ratio equal, then all entries of the /lg matrix, defined in 
(§, are real numbers. In addition if the Ybus matrix has diagonal elements that are non¬ 
negative real, non-positive imaginary and off diagonals are non-positive real, non-negative 
imaginary, then all entries of /lGj £^re non-negative real numbers. 

Proof. Let 
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, Yll = 

^nl 

bn2 

■ ^nn 
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and Ofej = Reiflkj) +ilm{akj), b^j = Re{b]^j) +ilm(bkj)- Then for each row k = 1 , 2 , ...,n 
of Ylg, Yll we have 


Im{akj) _ Imjbkj) _ 
Re{akj) Re{bkj) ^ 


k = 1 , 2 ,..., n. 


From (§ we have 

YllFlg = —Ylg, 

or, equivalently, 


■ &11 

bi 2 

bln 


Oil 

012 

dim 

^21 

^22 

^2n 

Flg = - 

021 

022 

0 ‘ 2 m 

bjii 

bn 2 

bnn 


Oyil 

an2 

■ ^nm 


or, equivalently. 


Re{bii)[l + iui] i?e(6i2)[l + iwi] 

Rclb2i)[l + iu2] i?e(&22) [1 + *'^ 2 ] 


Re{bin)[l + mi] 

-Re(&2n)[l +iU2] 


Flg = 


R^(,bnl^\^ F R^(,bn2^\^ F . . . Rc(^bjiji'j^l 


or, equivalently. 


Re{aii)[l + iui] 
Re{a2i)[l F iu2] 

Re{ani)[l F m„] 


Re{ax2)[l F mi] 
Re{a22)[l F m2] 

-Re(a„2)[l + iun] 


Re{aim)[l F mi] 
-Re(a2m)[l + iu2] 

Re{ 

^nm )[1 + m„] 


diag { 1 + mi, l + m2, 


1 + m„ } 


Re{bii) Re{bi2) 
i?e(&2i) ^e(&22) 


R^ib\ri) 

Re{b 2 n) 


Flg = 


R^ibnl') R^{bji2') ■ ■ ■ Rci^bjin) 


—diag{ iFiui, lFiu2, 


iFiUn } 


Re{aii) Re{ai2) 

Re{a2i) Re{a22) 


Rg{^Qiji2^ 


RG{cL2m) 

RGi^Qjim) 


or, equivalently, since det{diag { 1 + iui, 1 F m2, • ■ •, 1 + m„ 1 ) 7^0 


Re{YLL)FLG = —ReiYLG)- 


This is a liner system where the known matrices have only real entries. Hence, any solution 
for Flg lies inside of i.e. all entries of Flg are real numbers. In addition since 
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the Ybus matrix is assumed to have diagonal elements that are non-negative real & non¬ 
positive imaginary and off diagonals are non-positive real & non-negative imaginary, the 
matrix ReiX^o) has all its entries non-positive, i.e. the matrix Rc^lg) has all its entries 
non-negative and the matrix Rc^Yll) has its diagonal elements non-negative and its off 
diagonals non-positive, i.e. the inverse (or pseudo inverse) of Rc^ll) has all its entries 
non-negative, see [22] • Hence the matrix will have all its elements non-negative since 
it is given from the product of two matrices, the inverse (or pseudo inverse) of Re{Yi^]b) 
and the matrix Re{YLG) ■ The proof is completed. 


If we use a short line model of the system and so neglect shunt elements, then each 
row of the Y^us matrix sums to zero. If we include shunts, then each row sums close to 
zero. We can state the following Theorem. 


Theorem 2.1. Assume the Y^us matrix as defined in Q and the matrix F^q, as defined 
in Then the row sum of Flg is one if the rows of Y^us sum to zero and detXLh) ^ 0 
and is close to one if the rows of Ybus sum close to zero or sum to zero and det{YLi,) = 0. 


Proof. Let for system 0 


and 


Where 


Ylg 




yll = 


(7) 

( 8 ) 


Oil 

«12 

dim 
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We assume that each row of Y^us sums to zero, i.e. V row z = 1,2, ...,n of the matrices 
0, 0 we have “y + YTj=i equivalently 


From 0 we have 


'y ' — y ' — 1, 2,..., 7Z. 

i=i i=i 

YllFlg = —Ylg, 


(9) 

( 10 ) 


Let 


flg = (11) 

Then Vz = 1,2, ...,n by substituting 0, 0, ( [TT| ) into ( [Tq| ) and by calculating every 
component in row z, we get 


1 ^ik^kl 
— ^^k—1 ^'ikC^k2 


^im — ^^k—1 ^ikCkm- 
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By taking the sum of the above equalities, V* = 1, 2 ,n we arrive at 


n n n 

— [(lil + 0-i2 + ■■ + dim] = bikCk2 + ’ ' ’ + bikCkm, 

k=l k=l k=l 


or, equivalently, 

m m n 

- ^ikC^kj • 

j^l j^l k^l 

By using ([^ on the above expression we get 


X! - X! X! ^ik^kj j 

j = l k^l 

or, equivalently, 

n m n 

^ik^kj 1 

k—1 j—1k—1 

or, equivalently, 


m 

H“ ^i2 H“ ' ' ' H“ ^ H“ ^12^-2j H“ ' ' ' H“ 

i-1 


or equivalently, 


+ ^22 + * * * + bin — biiCii + bi2C21 + • • • + ^mCnl + 
^2lCl2 + bi2C22 + • • • + binCn2-\- 


bilCim + bi2C2m H“ ' ’ ’ H” bin^nm i 


or, equivalently. 


bil + bi 2 + ■ ■ ■ + bin — ^il(cil + Ci 2 + • • • + Cim) + 
bi2{C2l + C22 + • • • + C2m) + 


bini^djil Cn2 “t“ * * ’ “t“ Cnm') ^ 


or, equivalently, 

mm m 

bil “t“ bi2 “h * * * H“ bin — bii{ ^ bi2 ) H“ * ’ ‘ “h 

i=i i=i 


or, equivalently. 


mm m 

“ 1] + baiC^ C2j) - 1] H-h binlC^ Cnj) - 1] = 0. 

i=i i=i i=i 


( 12 ) 
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Let Vi = 1, 2,n 


Qi — ( '^ ' ^ij) 1 - 
i=i 


Then by replacing (13) into (12), Vi = 1, 2,n we have 


bilQl + bi2q2 + • • • + binQn — 0, 


(13) 


or, equivalently, 


[ bii bi2 • ■ ■ biji j 


91 

92 

Qn 


= 0, Vi = 1,2,..., n. 


Since the above expression holds Vi = 1, 2,..., n we have equivalently 


bii bi2 
^21 ^22 

^In 

^2n 


qi 

92 


1 - 

o o ■ 

1__ 

bnl bn2 

^nn 


Qn 


_ 0 _ 


By setting q — 

9i 

1 0n,l — 

■ 0 ■ 
0 

and by using ([^ we have 


Qn 


0 



^LL ' Q — 0n,l ■ 


From the above expression if detiY^i^) ^ 0, then q = 0n,i, or, equivalently, qi = 0, 
Vi = 1, 2,..., n, i.e. by using ( [T^ 

ra 

E"V = 1- (14) 

i=i 


If detiY^i) = 0, then q = 0^^ (by using the pseudo inverse of Fll via the SVD method), 
or, equivalently, qt = 0, Vi = 1,2, ...,n, i.e. by using (13) 





(15) 


Thus, from (14) and (15) every row of Fi^q will sum to 1 if det{Y]^i) ^ 0 and will sum 
closely to 1 if det{YLL) — 0. If the rows of Ybus sum close to zero, then with similar steps 
we arrive at The proof is completed. 
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Conclusions 


The partitioning of the matrix has opened numerous fruitful avenues in power system 
analysis, and it is hoped that a clearer understanding of the matrix properties underpinning 
these partitions will support this strand of research. This work has proved two numerically 
observed properties of the Tlg matrix: future work may consider the matrix characteristics 
of the other sub-matrices derived by the Ybus partitioning approach. 
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